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We present the latest results on the study of atmospheric neutrino oscillations with the MACRO detector at 
Gran Sasso. Two sub-samples of events have been analysed, both in terms of absolute flux and zenith angle 
distribution: high energy events (with {E^) ~ 50 GeV) and low energy events (with {Ev) ~ 4 GeV). The high 
energy sample has been used also to check the ^sterile oscillation hypothesis and to estimate neutrino 

energies using Multiple Coulomb Scattering (MCS) informations. All these analyses are mutually consistent and 
strongly favour the <-> Vt oscillation hypothesis with maximal mixing and Am^ = 2.5 • 10“®eV^. 


1. INTRODUCTION 

The MACRO detector 0 was located in the 
Gran Sasso Laboratory (Italy) and allowed the 
study of atmospheric neutrinos detecting upgo- 
ing muons produced in CC interactions inside or 
around the detector. It used scintillation coun¬ 
ters for tagging events by time-of-flight (TOP) 
measurements, and limited streamer tubes for 
tracking. Three categories of z/^-induced muon 
events were studied (see Fig. 1) upthrough- 
going muons, produced in the rock below the de¬ 
tector and crossing the whole apparatus, with 
an average energy {E,j) ~ 50 GeV; 2) internal 
upgoing (lU) muons, produced inside and leav¬ 
ing the detector from above; 3) upgoing stopping 
muons (UGS), produced below and stopping in¬ 
side the detector plus internal downgoing (ID) 
muons, produced inside and leaving the detector 
from below. Since at least two scintillation coun¬ 
ters are needed to perform a TOF measurement, 
the events of samples UGS and ID are studied 
on a topological basis and are indistinguishable. 
Therefore, they are studied together. The average 
neutrino energy for samples 2 and 3 is about ~ 4 
GeV, and we collectively refer to these events as 
to the Low Energy (LE) sample ||] , while events 
from sample I belong to the High Energy (HE) 
sample |Q. 


*List of authors and institutions: see Ref ll 


2. HIGH ENERGY SAMPLE 

Events in the HE sample have been selected by 
measuring the TOF between two layers of scin¬ 
tillators. For events crossing three layers of scin¬ 
tillators, a linear ht of the times as a function 
of the path lengths has been performed, in or¬ 
der to reduce possible fake events. A detailed 
tuning of the TOF measurements provided a re¬ 
duction factor of ~ 10^ of the downgoing atmo¬ 
spheric muons. Other sources of background have 
been reduced by applying a cut on the match¬ 
ing positions between scintillators and streamer 
tubes. Soft upgoing hadrons coming from pho- 
tonuclear interactions of muons outside the detec¬ 
tor could mimic an upgoing muon: to reduce this 
background we require that each upgoing muon 
crosses at least 200 g/cm^ of absorber in the lower 
part of the detector Q . Muons coming from par¬ 
ticular azimuthal regions, where the rock amount 
is not sufficient to reduce the number of down¬ 
going muons, have been discarded. 863 events 
survived these cuts (809 after background sub¬ 
traction), corresponding to a livetime of 6.16 yrs. 

Monte Carlo (MC) simulation of upgoing muon 
events has been performed using the v flux com¬ 
puted by the Bartol group Q , neutrino cross sec¬ 
tions from Ref. and muon energy loss in the 
rock from calculations of Ref. [Q. The overall 
theoretical uncertainty on the absolute upgoing 
muon flux is 17%, while the systematic error on 
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Figure 1. Schematic picture of the different 
topologies of r^-induced events in MACRO. From 
the left: 1) upward throughgoing muons, 2) in¬ 
ternal upgoing (lU) muons, 3) upgoing stopping 
muons (UGS) and 4) internal downgoing (ID) 
muons. Samples (3-1-4) are studied together, see 
text. 



Figure 2. Angular distribution of upthroughgoing 
muons. Black points are data, with statistical and 
systematic errors summed in quadrature. Shaded 
region is Monte Carlo prediction, assuming no os¬ 
cillations and with a 17% error on the scale. The 
solid line is the Monte Carlo prediction assuming 
t'/p i't oscillations with maximal mixing and 
Am^ = 2.5 ■ IQ-^eV^. 


the shape of the angular distribution is 5%. 

Compared to MC predictions the number of 
measured events is small and the shape of the 
angular distribution is different, Fig. |[ The ra¬ 
tio between the observed number of events and 
the MC prediction is Rhe = 0.721 ± 0.026stat ± 
0.043sys ± 0.123tft,. The shapes of the experi¬ 
mental distribution and MC predictions assum¬ 
ing no oscillations, analysed in terms of give 
an agreement probability of 0.2%. Assuming 
^ Vt oscillations, the angular distribution and 
the absolute number of events give a best fit for 
maximal mixing and Am^ = 2.5 • 10“^eV^, with 
a probability of 66%. The result of the fit is 
shown as a solid line in Fig. H. 

2.1. <-> Vsterile oscillations 

The weak potential of Vf. and Vsteriie with mat¬ 
ter is different from the and Vr potentials. 
On the other hand, and Vt potentials are 


equal: this difference translates in a distortion 
of the oscillation pattern (and hence of the an¬ 
gular distribution) in the ^ Vsteriie oscillation 
with respect to the <-> Vr case. Matter ef¬ 
fects become important when E^/\/S.w?\ > 10“^ 
GeV/eV^, i.e. for HE events. A detailed analy¬ 
sis [|| has shown that the best estimator to dis¬ 
entangle the two hypotheses is the ratio R be¬ 
tween the number of events with -1 < 0 < - 
0.7 and the number of events with -0.4 < 0 < 
0, where 0 is the zenith angle. In this ratio, 
most of the theoretical uncertainties cancel and 
the overall uncertainty, combined with the exper¬ 
imental systematic error, is ~ 7%. The measured 
ratio is i? = 1.48 ± 0.13stat ± O.lOsys- This result 
is shown in Fig. ^ together with the MC pre¬ 
diction, as a function of the Am^. In the best 
fit point, the expected values for z/^ ^ Vt and 
Vsteriie are Rt = 1.72 and R^- = 2.16 re¬ 
spectively. The ratio of the probabilities to obtain 
values as low as the observed one is Rprob = 157, 
therefore Vp Vsteriie can be excluded with re- 
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Figure 3. ratio of the number of events with -1 
< 9 < -0.7 with respect to the number of events 
with -0.4 < 0 < 0. The black point is the mea¬ 
sured value (with error bar), the solid lines are the 
prediction for <-> Vr and i^steriie cases. 


spect ^sterile at 99% CL. 

2.2. Energy estimate by means MCS 

The limited streamer tube system of the 
MACRO detector has been used in order to es¬ 
timate muon energies (and hence neutrino ener¬ 
gies) using MCS informations. A first analysis, 
which used streamer tubes in digital mode (with 
a spatial resolution of cr ~ 1 cm), successfully 
showed the feasibility of the method |^. In or¬ 
der to improve the sensitivity of the analysis to 
higher neutrino energies, we used the streamer 
tubes in drift mode, improving the resolution by 
a factor ~ 3.5 0. Two dedicated test beams at 
CERN PS/SPS checked electronics, detector per¬ 
formances and analysis tools. Starting from MCS 
sensitive variables, a MC trained neural network 
was used for muon (and neutrino) energy recon¬ 
struction. These results allowed to increase the 
sensitivity of the analysis up to neutrino ener¬ 
gies of ~ 100 GeV. We separated the HE sample 
in four different energy regions: the results show 
that the agreement between data and MC pre¬ 
dictions follows the energy dependence expected 
in the oscillation hypothesis with the parameters 



Figure 4. Ratio DATA/MC{no — osc) as 
a function of the estimated logioL/E„ (black 
points). The solid line is the ratio MC{siri^26 = 
1, Am^ = 2.5- 10“^eV^)/MC'(no —osc), with cor¬ 
responding systematic error (see text). The black 
square refers to lU events. 


given in the previous section ]ll|] . 

The distance travelled by neutrinos inside the 
Earth (reconstructed with a precision of ~ 3%) 
was used to estimate the L/Ei, on a event-by- 
event basis. This is shown in Fig. where the 
ratio DAT A/MC {noose) as a function of the es¬ 
timated logipL/Ei/ is plotted. Errors include un¬ 
certainties both on the flux and on the shape. 
MC prediction is also reported, with the same pa¬ 
rameters of the standard analysis, showing that 
the trend of the data is the one we expect in the 
Vfi ^ Vt oscillation hypothesis. 

3. LOW ENERGY SAMPLE 

lU and ID-I-UGS events are produced in CC 
interactions inside the detector, with a small con¬ 
tamination of NC and z/g events (^ 13% and ^ 
10% respectively). lU events have been selected 
by TOF measurement and by topological criteria, 
e.g. the requirement of the interaction vertex in 
the fiducial volume of the detector. The ID-I-UGS 
sample has been selected only with topological 
criteria. Since no timing informations were avail¬ 
able for these events, particular attention has 
been devoted in order to reject events not per¬ 
fectly confined in the detector. These goals have 
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Figure 5. Angular distribution of lU events (a) 
and ID+UGD events (b). Black points are data, 
with statistical and systematic errors summed 
in quadrature. Dashed lines are Monte Carlo 
predictions, assuming no oscillations and with a 
17% error on the scale (shaded region). Solid 
lines are the Monte Carlo predictions assuming 
Vt oscillations with maximal mixing and 
= 2.5 ■ IQ-^eV^. 


been reached also by means of visual scanning 
with the Event Display. After these selections, we 
remained with 154 and 262 events for the lU and 
ID+UCS, corresponding to 5.8 yrs and 5.6 yrs 
of livetime respectively. In Fig. || we show the 
angular distributions of these two samples and 
the MC predictions with and without the oscil¬ 
lation hypotheses. A conservative value (25%) of 
the theoretical uncertainty on the absolute scale 
is reported. There is a good agreement between 
data and MC prediction with the oscillation pa¬ 
rameters given by the HE sample analysis. 

4. CONCLUSIONS 

We presented different and independent anal¬ 
yses performed on muon events induced by at¬ 
mospheric neutrinos. All the analyses are mutu¬ 
ally consistent and strongly favour the 
oscillation hypothesis with maximal mixing and 
Am^ = 2.5 • 10“^eV^. Absolute flux measure¬ 
ments, shape distributions and energy estimate 
by MCS were used in order to compute allowed re¬ 
gions in the space of oscillation parameters. This 
is shown in Fig. ^ where the corresponding con¬ 
tours are computed according to a analysis 
using the prescriptions of Ref. |Q. A global 
analysis of all the samples is in progress. 
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Figure 6. 90% CL contours for allowed region 
in the {sin^20, plane. Curves 1 and 2 refer 

to upthroughgoing muons using flux-|-angular dis¬ 
tribution and energy estimate respectively. The 
curve 3 refers to the LE sample. 
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